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NOMENCLATURE 


CA] 

AIC 

{a} 

c 

g 

CSA3 

Ci] 

i 

m 

m 

k 

CM] 

m 

{U} 

V 

X 


aerodynamic influence coefficient matrix 
abbreviation for aerodynamic influence coefficients 
vector of modal participation factors 
reference chord length 

scalar approximation of the structural damping of the system 

matrix of generalized aerodynamic forces 

identity matrix 

imaginary constant 

stiffness matrix 

diagonal generalized stiffness matrix 
reduced frequency parameter = (oc/2V 
inass matrix 

diagonal generalized mass matrix 
structural displacement vector 
velocity 

eigenvalue of the complex eigenmatrix ( = w /(1+ig)) 
matrix of normal modes from the vibration analysis 
circular frequency 


Subscripts 

r refers to rigid body modes from the vibration analysis 

f refers to flexible modes from the vibration analysis 

i ith component of a vector or the ith column of a matrix 

n. Superscripts 

( matrix inverse 

Note: Additional notation is defined in the text. 
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FLIJT - A PROGRAM FOR AEROELASTIC STABILITY ANALYSIS 


Erwin H. Johnson* 

Ames Research Center 

SUMMARY 

FLUT is a computer program that can be used to evaluate the aeroelastic 
stability of aircraft structures in subsonic flow. The algorithm synthesizes 
data from a structural vibration analysis with an unsteady aerodynamics 
analysis and then performs a complex eigenvalue analysis to assess the system 
stability. This document is divided into two main parts. The first of these 
describes the theoretical basis of the program. Special emphasis in this 
section is placed on some innovative techniques which improve the efficiency 
of the analysis. 

The second section provides the user information needed to efficiently 
and successfully utilize the program. In addition to identifying the required 
input, this section summarizes the flow of the program execution and points 
out some possible sources of difficulty. The use of the program is demon- 
strated with a listing of the input and output for a simple example. 


*National Research Council Postdoctoral Research Associate. 


I. Introduction 


The Advanced Vehicle Concept Branch at the NASA-Ames Research Center has 
been charged with responsibility of assessing new aircraft concepts and pro- 
viding quantitative data on how these concepts might best be implemented. 

Many of these concepts were characterized by lightweight high-strength materials 
and high aspect ratio configurations. This naturally led to concern about the 
aeroelastic stability properties of the concepts and a need was recognized for 
a technique that would perform this stability analysis. In particular, there 
was great uncertainty as to the aeroelastic stability properties of a subsonic 
oblique wing transport concept under study by the Locicheed-Georgia Company. 

Linder these circumstances, ahe FLUT program was developed with a special 
emphasis placed on using an aerodynamics package that did not require aircraft 
symmetry about the fuselage. It is felt that this program can be useful in a 
number of applications, which therefore motivates this documentation. 

The first part of this document describes the analysis on v/hich the program 
is based. Special emphasis is given to features of the program that may be new 
to workers in the field. It does not attempt to provide an in-depth description 
of flutter analysis techniques; a more complete description of such techniques 
can be found in any of several texts on aeroelasticity (e.g.. Refs, i and 2 J. 

The remainder of this report is divided into three sections. The first 
defines the input required by FLUT, The next provides some information on 
program execution, and describes its operation with the use of a flow chart. 

The last section presents a simple example. 

Additional information on the program usage is contained in the source 
listing. The majority of the variables used are explicitly identified in 
comment cards and the program flow is outlined. The unsteady aerodynamics 
package contained in FLUT was developed by Giesing, Kalman and Rodden and 
is explained in detail in References 3 and 4 . 

It must be noted that much of the basis for the methods of analysis used 
in this program were obtained from a program by Crittenden and Weisshaar, 

(Ref.. 5 ). 
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II. Basic Matrix Representation 


The program FLUT is based on what has come to be known as the k method 
of flutter analysis (Ref. 1) which starts from the matrix equation of motion: 

[0+ig)[K] - + [A])] {u} = 0 CD 

This equation has degrees of freedom that are che actual structural displace- 

^ • 

ments. The mass and stiffness representations are therefore the finite element 
representations developed by any computer code similar to one described in 
Reference 6 „ 

The aerodynamic matrix contains the coupling between structural displace- 
ments and the aerodynamic forces they create. In FLUT, ^nis matrix is calculated 
using doublet lattice methods contained in routines documented in References 3 
and 4. This matrix is a function of the free stream Mach number and the reduced 
frequency parameter k. The point to be stressed here is that aerodynamic 
matrices of the form required by Eq. (1) are actually calculated during 
program execution. This contrasts with other programs which calculate 
generalized aerodynamic forces directly. More v/ill be said on this in the 
"Efficiency Concepts" section. 

The matrix eigenvalue problem given by Eq. (1) has a large number of 
degrees of freedom and is not in a form needed for a complex eigenvalue solution. 
A more convenient form is one that uses normal modes to represent the displace- 
ments. This IS done by a transformation of the form 

mn 

Where mn is the number of normal modes retained in the analysis. If this 
transformation is placed in Eq, (1) and the result is premultiplied by 
the transpose of the matrix of the normal modes, the resulting equation is of 
the form 
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Because normal modes were used, the generalized mass and stiffness matrices 
are diagonal. The generalized aerodynamics matrix is complex and nonhermitian. 

The complex eigenvalues of the matrix 

[XrkET/?] + [SA]) w 

are then the desired results. 

When rigid body modes of the aircraft are included in the analysis, the 
generalized stiffness matrix has zeroes on the diagonal. This makes it impossible 
to perform the inversion called for by Eq. [4). The next subsection pro- 
vides a means of dealing with this situation. 

Reduction of Rigid Body Modes .- When rigid body modes are present, 


Eq. (3) can be subdivided into the form: 



A reduced problem is obtained by solving Eq. (6) for and placing the 
result in Eq, (7). 


This equation is now well behaved and its solution can be obtained in the 
same way as one without rigid body modes. Although this technique for reducing 
rigid body modes is not generally known, it has been previously discussed by 
Cross and Albano in Reference 7. 

Equation (4) or the system resulting from Eqs. (6) and (7) can now be evalu- 
ated to find the complex eigenvalues X and their corresponding eigenvectors 
(which represent modal participation factors). The complex eigenvalue analysis is 
performed using routines extracted from Reference 8 which includes descriptions of 
the algorithms. Briefly, the routines transform the matrix to upper Hessenberg 
form and then use a "modified LR method" to extract the eigenvalues and eigenvectors 
With the solution obtained, the only remaining task is to get the results 

in a useful form for output. The eigenvalues contain the needed stability 

2 * 

information. The X's obtained represent w /(1+ig); therefore, the velocity, 
damping, and frequency information can be obtained from 

g = -Imx/Rex 

= (l+g^)ReX (8) 

V - 

Section VI provides added detail on how these results are evaluated. 

The eigenvectors from the eigenvalue analysis are the modal participation 
factors for the flexible models. It is possible to obtain the participation 
factors for the reduced rigid body modes through a relationship based on 

Eq* (6): 

GA^^]’’[GAf^] (a^) (9) 

If the actual structural deflections are required, they can, of course, be 

obtained from the superposition relationship of Eq. (2). 

The program actually solves for 1/X and derives g, w, and V from this 
variable. 


5 


Ill Efficiency Concepts 

The utility of any analysis can be enhanced by making it as efficient as 
possible. This is particularly true for a flutter analysis; since, typically, many 
repetitive analyses are required. Steps have been taken in FLUT to perform 
these repetitive analyses with a minimum of additional effort. Two steps that 
are of sufficient general interest to be described here are (1) aerodynamic 
interpolation and (2) optimal use of aerodynamic influence coefficients. 

Both of these concepts derive their importance from the fact that the 
time required to perform an unsteady aerodynamics analysis at a single set of 
k and M values can be as much as an order of magnitude higher than the time 
required to perform the complex eigenvalue analysis. It is therefore very 
important to keep the aerodynamics calculations to a minimum. 

Interpolation of the aerodynamics can be done to give high quality results 
with a minimum of effort, thereby minimizing the aerodynamic calculations. 

The interpolation is performed by fitting a quadratic in 1/k through 
generalized aerodynamics forces that were in turn calculated using aerodynamic 
influence coefficients (AIC's) computed at three reduced frequency values. 

The preceding sentence brings up several points: 

Cl) The AIC's are calculated using the routines of Reference 4. 

(2) Generalized aerodynamics, rather than AIC's were used in the fit 
since the matrices involved are smaller and there is no loss in 
accuracy, 

(3) The inverse of the reduced frequency was used to obtain the fit since 
it provided results vastly superior to those obtained using the 
reduced frequency directly. This is borne out by Table 1 which com- 
pares results for one branch of a V-g analysis obtained usjng either 
“exact" aerodynamics (aerodynamics computed directly without inter- 
polation), aerodynamics interpolated using 1/k as the independent 
parameter and those obtained using k. It can be seen that the results 
obtained using 1/k are within one percent of the exact values, while 
the results using k are entirely unacceptable. 
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Table 1.— Comparison of Aerodynamic Interpolation lechniques 



“Exact" 

Interpolation with j< 

Interpolation with 1/k 


V 

g 

(0 

V 

g 

0) 

V 

g 


.01* 

678.5 

.053 

.6785 

678.5 

.0526 

.6785 

678.5 

.0526 

.6785 

.03 

702.0 

-.113 

2,11 

336.1 

.073 

1.01 

702.4 

-.112 

2.11 

.05* 

718.3 

CM 

• 

1 

3.59 

718.3 

-.247 

3.59 

718.3 

-.247 

3.59 

.10 

720.4 

-.378 

7.20 

141.0 

.759 

1.41 

719.2 

-.379 

7.19 

.15 

656.8 

-.380 

9.85 

86.4 

.628 

1.29 

655.8 

-.382 

9.84 

.20* 

573.2 

-.348 

11.50 

573.1 

-.348 

11.46 

573.2 

-.348 

11.46 

.30 

426.2 

-.278 

12.79 

12.53 

+ .02 

.375 

428.2 

-.269 

12.85 

.50 

263.1 

-.175 

13.16 

3.18 

-.025 

.159 

265.7 

-.152 

13.29 

, 


* The aerodynamics used for the interpolation were computed at these points. 
Therefore, the interpolated results at these reduced frequencies are the 
exact results. 











(4) Presently only interpolation in 1/k is performed. It may be desirable 
to include interpolation in Mach number for some applications. 


Finally, it should be stressed that the interpolation factors that weight 

the contributions from each of the three basic general izeH aerodynamic matrices 

2 

are the same for all the mn elements of the aerodynamics matrix, i.e., 

[GA^nt^ = fiCGA^] + f2 [GA2] + f3 CGA3] ( 10 ) 


Here, int refers to interpolation results, f -j , f 2 , and f 3 are constant scalars 
which weight the three exact generalized aerodynamics matrices. 

I - (lu-b-vlttKylti-vka) 



(vlti - vtsKvl 

li-vlti) 

(ulti- 

It; - IrlfLz) 



where the convention v 


k. 

J 


= 1/kj is used. 


( 11 ) 


The second efficiency, the use of AIC's, allows a number of different 
vibration analyses to use the same aerodynamics results. This contrasts with 
unsteady aerodynamics programs that compute the generalized aerodynamic forces 
directly which require the calculation of new forces each time a different vibration 
analysis is made. Occasions where the AIC's can save computer resources are: 

(1) When structural changes in the aircraft are being studied. 

(2) When various boundary conditions are being evaluated, e.g., it is 
sometimes desirable to see what the effect of rigid body roll is 
on the flutter analysis. The AIC's remain the same while the mode 
shapes differ drastically from a clamped wing condition. 
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(3) When investigations into the sensitivity of the flutter* speed to 
parameters such as static unbalance, attached masses* and fuselage 
inertias are being performed. 

IV Program Input 

The following is a list of inputs and their associated format statements. 
Where additional clarification is required, reference is made to notes con- 
tained at the end of this section. The input required for the aerodynamics 
package (Blocks C through 0)was set by the algorithm's developers and is 
explained in Reference 4. 

Dimensions used in the program are arbitrary, with the exception that they 
must be self-consistent; e,g», deflections for the mode shapes cannot be in 
inches while the chord length is in feet. 

The format of this description of data is the following: A letter signi- 
fying the data block is given in the left hand column. The input is then 
listed and Its associated format is given in parentheses in the right-hand 
column. A description of the data is then given before the next data Mock 
is listed. 

A. TITLE (20A4) 

Title is any eighty character title. 

B. N',HN, NA, NP, NAR, NRD, NRFQ, NTRF, NREAD, NW, IPRl, 1PR2, NAERO (2013) 

N - Total number of unconstrained iegrees of freedom from 
the vibration analysis (N<90) 

MN - Number of normal modes retained from the vibration 
analysis (MN<20) 

NA - Number of aerodynamic degrees of freedom (NA<100) 

MR - Number of rigid body modes included in the vibration and 

flutter analysis (NR<3) 

NAR - Number of structural degrees of freedom that are not at 

an aerodynamic grid point and must therefore be eliminated. 

See Note a. (NAR<10) 

NRD - Should be zero. See program listing. 

NRFQ - Number of reduced frequencies for which the aerodynamics 
are calculated explicitly (NRFQ<3) 

NTRF - Number of reduced frequencies at which the complex eigen- 
value problem is to be solved. See Note b. (NTRF<20). 
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NREAD 


Number of the tape on which the vibration analysis results are 


NW 


IPRl 


IPR2 


NAERO 


written. This number must, of course, be the same as the tape 
number attached in the JCL. 

A tape number for writing plots data. It is recommended that 
NW = 0 be used so that no plots data will be generated. The 
plots information is unsatisfactory. 

Print command. If it is greater than zero, the vibration analysis 
input, the generalized aerodynamics, the complex eigenmatrix and 
the modal participation factors are printed. iPRl = 0 
suppresses these prints. These are intermediate prints, intended 
for use during program debug. 

- A second print command. If it is greater than zero, the aerodynamic 

influence coefficients and the flutter mode shapes are printed. 

IPR2 = 0 suppresses these prints. The outputs from this 
parameter can be considerable and it is recommended that it 
usually be set to zero. 

- Parameter indicating how the aerodynamics are to enter the 

formulation. 

= 1 - Program computes the AIC's at the NRFQ reduced frequencies 
and saves results on tape 11. (Of course, the 
results are actually only saved when the appropriate 
cataloging card is included in the JCL.) 

- 2 - Previously computed aerodynamics are read from tape 11. 
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If (NAERO = 2) go to Block P. 

C. TITLE 

This is a second title used for the aerodynamics package. 

D. FMACH, ACAP, FL. B2, NDELT, NP. NB, NRF, NCORE, Nl, N2. N3, M4 (4F10.0, 

412, 16, 411). 

FMACH - Free stream Mach number of the aerodynamics calculations. Must be 
the same as AMACH in Block R. 

ACAP - Reference area, total area of the wing 

FL - Reference chord. An average chord of the wing. Must agree with CBAR 
in Block R. 

B2 - Wing semispan 

NDELT - Symmetry flag for the ZY plane 
= 1 Symmetric aerodynamics 
== -1 Antisymmetric aerodynamics 

=0 No symmetry (this option is used for obliquely configured wings, 
for example.) 

NP - Total number of aerodynamic panels. See Note c for a definition of the 
terms panel, strip and box. 

NB - Total number of bodies. FLUT has never been used with bodies and would 
require some reprogramming to do so. Therefore, NB should be zero. 

NRF - Number of reduced frequencies for which the aerodynamic analysis is 
to be performed. NRF must agree with the value of NRFQ in Block B. 

NCORE - This parameter is a flag to determine which set of routines are to 
be used for solving linear systems of equations. See Itite d. It 
should be set equal to the total number of boxes times the value 
of NA from Block B. 
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N1 - A flag which should alv/ays be zero for this configuration 

of the program. 

M2 - A flag delineating v;hether polynomial modes on AIC's are 

being calculated. Should always be 1 for this configur- 
ation. 

M3 - If N3 is >0, voluminous intermediate results are printed 

out. It is recommended N3=0 be used unless detailed debug 
is being done, 

N4 “ Print parameter for slender body data. N4 should always 

be zero for this configuration of the program. 


E. (RFREQ(I). 1=1, NRF) _ (6F10.0) 

filC 

RFREQ - Vector of reduced frequency values = k = . These 

numbers must be the same, and in the same order, as those 
read in by Block P. 

Repeat input Blocks F through I, for each panel. 


F. XCAP(l), XCAP(2), XCAP(3), XCAP(4), YCAP(l), YCAP(2) (6F10.0). 

Blocks F and 6 read in panel coordinates as shown in Figure 1. The x 
coordinate is in the streamwise direction while the y coordinate is per- 
pendicular to the flow. The streamvnse edges of the panels are parallel 
to the flow. The origin of the coordinate system is arbitrary. 

YCAPO)< YCAPC2) 

XCAP(1)< XCAP(2) 

XCAP(3)< XCAP(4) 

G. ZCAP{1), ZCAP(2), NS, NC, COEFF, (2F10.0, IX, 213, 3X, FIO.O). 

ZCAP - Vertical locations of the streamwise panel edges. See 

Figure 1. 

NS - Number of spanwise division boundaries for the panel. 

It is one greater than the number of strips. 

NC - Number of chordwise division boundaries for the panel. 

It is one greater than the number of boxes per strip. 
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X, 


COEFF - Scale factor for body deflections. It should be set to 1. 

H. (TH{I), 1=1, KC) {6F10.0) 

TH{I) = - Chordwise divisions of the panel in fractions of the chord. 

Usually, e-jis 0-0 and is 1.0. 

I. (TAU(I), I = 1, NS) (6F10.0) 

TAU(I) = T^. - Spanwise divisions of the panel in fractions of the panel span. 

Usually, T is 0.0 at the left edge of the panel and is 
1.0 at the right edge. 


J. NSTRIP, NBl, NB2, NPRl, NPR2, NPR3, JOBNO, J0BN02, NGUST, OSPECS, NPC, NSV, 
NBV, NYAW (IX, 313, IX, 313, 115,. 15, IX., 513, 212) 


NSTRIP 

NBl 


NB2 


NPRl 


NPR2 


NPR3 


Total number of strips. See Note c. 

Sequence number of the first strip. For tnis configurations 
strip numbering starts from the left wing tip (or the inboard 
edge if only the right wing is being considered) 

Sequence number of the last strip. For this configuration, strip 
numbering ends at the right wing tip. 

Print flag. If NPRl = 1, intermediate results are presented. 

It is recommended that this be zero so that printing is suppressed. 
Print flag. If NPR2 = 1, AlC's are printed. During early runs, 
these can be interesting, but it is recommended the NPR2 = 0 
be used to suppress print and limit printout. 

Print flag for stability derivatives 
= 0-No stability derivative information is calculated 
= 1 - dynamic stability derivatives are printed 
= 3 - static and dynamic derivatives are printed. For this 
case, NRF must be > 2 and RFREQ(l) must be 0.0. 
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It is recommended the NPR3 1 be used during early runs of 
a study since the derivative information provides a qualitative 
check of the program. During production runs, NPR3=0 should be used. 
NPR3 =3 cannot be used in a flutter study since the AIC’s at 
RFREQ - 0 are not compatible with the aeroelastic stability 
algorithm. See rtite e. 

JOBNO " Tape number for AIC’s. Use any four digit number. It appears 
that this input is never used. 

J0BM02, riGUST, dSPECS, NSV, NBV, and NYAW are not used for this configuration. 

Therefore their fields can be zero or blank. See Reference 2 for 
their definition and possible application. 

NPC - Mode selector for AIC generation. ' °C=0 gives plunqe, pitch, 
control surface and tab degrees of freedom. NPC=1 gives 
“cambering" modes plus control surface and tab degrees of 
freedom. MPC = 0 should always be used for this configuration. 

K. (LIH(I, 1), LIMCI.2), LIM(I,3), I « 1, NSTRIP) (6{IX, 313)) 

LIMCI, 1) “ Number of the first box for strip I. 

LIM(I, 2) - Number of the last box for strip I. (See Note c.) 

LIM(I,3) - Probably for control surfaces. In any case, it is to be 
left blank in this configuration. 

L. (NOP(I), (IS(I, d), d = 1, 7), I = I, NP) (6(1X, 12, 711)) 

NOP(I) - Number of the panel 

IS (I, d) - Degree of freedom selector for panel. For this configur- 
ation, set IS{I, 1) = 1 for plunging. IS (I, 2) = 2 
for pitching and leave the rest blank. 

Repeat Blocks M and N for each panel. 


U 


M. 


XHI(I, 1) 


(FIO.O) 



i 


Tfi6 location of the reference axis (the line about which the aero- 
dynamics are collocated) at the left edge of the panel as a fraction of 
the chord. 

N. XH0(I, 1) (FIO.O) 

The location of the reference axis at the right edge of the panel 
as a fraction of the chord. NOTE: Other axes that could be read in by 
this card but are not relevant to this configuration, are not described 
here. 

O. NMD, NTA, N5, N6, N7, NTPl , NTPM, NMTP, NMTB (512, 2110, 212, 213) 
NTA, N6, NTPM, NMTP and NMTB can be ignored for this configuration, i.e., 
their fields should be left blank. 

NMD - Number of aerodynamic modes. This must agree with NA from 
Block B. This number is required only for NCORE ^ 3700. 

N5 - Flag to save intermediate results. Can be always set to zero 
for this configuration. 

N7 - Flag for pressure forces and moments. 

=1 calculate and print pressures 
=0 suppress calculation 

It is recommended M7 = 1 be used for early runs as a quali- 
tative check of 0 itput, then set to zero to limit output. 

NTPl - When NCORE is > 3700, an arbitrary number (say 720) must be 
used here to save intermediate results. Otherwise, the field 
can be left blank. 

This is the end of the aerodynamics package input. 

P. (RFQ(I), I == 1, NRFQ) (6E12.4) 

Values of the reduced frequencies at which the aerodynamics have 
been computed. These values must agree with the data in Block E. 
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s 


Q (RFT(I), 1 = 1, MTRF) 


(6E12.4) 


Values of the reduced frequency at which the complex eigen- 
value problem is to be solved. See N^te b. 

R. NRHO, Af^CH, GEAR, UCl , UC2 (13, 4E12.4) 

NRHO - Number of different densities at which the complex eigenvalue 

problem is to be solyecj. 

- Mach number at which the aerodynamics are calculated. This 
must agree with FMACH read in by Block D. 

GEAR - Reference chord length. Must agree with FL in Block D. 

UCl - Units conversion factor. UGl is the factor required to change 

the input length dimension to feet (e.g., if input is in 
inches, UGl = .083333) 

UG2 - Units conversion factor UC2 is a factor required to change 
the mass dimension from slugs to the units of the input. The 
only situation where this should differ from one is if kilograms 
are used, in which case UG2 = 14.594. 

S. IF(NAR.NE.O) READ (JRF(I) ,I = 1 ,NAR) (1013) 

JRF(I) = The integer number of the i degree of freedom that is 

not at an aerodynamic grid pointy See Note a. JRF(I+1 ) > JRF(I) . 

T. (RHONd), H(I), 1=1, NRHO) (2E12.4) 

RH0N{I) = p/po - Ratio of the i atmospheric density at the altitude 

the flutter analysis is to be performed to the sea level 
density. 

4* Vi 

H{I) - Altitude corresponding to the i^^ density (currently, this 

dimension is only in feet). This number is only used for output. 
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This completes the description of the data that is input from cards. Data 
from the vibration analysis and from the unsteady aerodynamic analysis is read 
from tapes. The aerodynamic calculations are all done within the program and 
are therefore self-consistent. The vibration results j, however, must be input 
in a very specific way: 

The vibration data is input in two steps, first the eigenvalues are 
read in and then the eigenvectcYS . The statements used for these reads are: 

Al. (LBDA(I), 1=1, MN) (unformatted) 

LBDA(I) - ith eigenvalue of the vibration analysis. This is the square 

2 

of the ith natural frequency. Units are (rad/sec) . 

LBDA(I + 1) > LBDA(I) 

Bl. ((PH(I,J), J = MN), I = 1, N) (unformatted) 

f'h 

PH(I,J) - The ith component of the j normal mode. 

Note, a. contains information on what the structural components are and 

how they relate to the aerodynamic degrees of freedom as v^iell as how to 
modify the program to accept other forms of the eigenvalues and eigenvectors. 

V. Notes 

These notes are intended to clarify and amplify on the descriptions of the 
input givm above. Most of the notes describe quirks of the program that could 
prove troublesome to a new user. 

(a) Structural and aerodynamic degrees of freedom. - 

Figure 2 shows a representation of an oblique wing planform with a set 
of grid points and aerodynamic boxes. Points 1 tarough 9 are structural 
grid points while points 2, 3, 4, 6, 7 and 8 are aerodynamic grid points 
as well. When a vibration analysis is being run to determine the natural 
frequencies and the normal modes, consideration must be taken of the 
aerodynamic representation that is to be used so that the requisite 
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degrees of freedom are obtained. This is because structural nodes must 
be placed at the centers of the aerodynamic boxes. Each of these grid 
points must include three, and only three, degrees of freedom: (1) 


deflection in the z direction, * (2) rotation about the global 
X axis, 9 , and (3) rotation about the global y axis, e . The dof's 
must be in the order given in the preceding sentence and must start from 
the left grid point and proceed to the right until the last grid point 
is reached* The eigenvectors must be normal izedv to give a unit mass 
matrix. 

Assume in Figure 2 that the structural model is represented by the 
nine nodes shown in the figure connected by eight beam elements. Each 
of the nodes has the three dof's listed above except node 5 which, in 
this example, has freedom only in the Z direction (allowing for rigid 
body plunge. The structural dof's are then: 

Structural dof Node Direction 



11 4 

12 4 

13 5 
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Structural dof 

Node 

Direction 

14 

6 


• 

* 

« 

• 

* 

• 

22 

8 

®y 

23 

9 


24 

9 


25 

9 

i. 

y 


There are, therefore, 25 structural degrees of freedom whose de- 
flections must be passed to FLUT in the order listed above. Of these 
dof's, several (numbers 1, 2, 3, 13, 23, 24, and 25) are not at aero- 
dynamic grid points. These are the JRF's given in Block S. 

As a further note on degrees of freedom, the aerodynamics package 
calculates forces due to pitch and plunge at the center of each strip. 
No forces are calculated due to 0 deflections and the plunging degree 

A 

of freedom is in the -z direction. Internal to the program, the 0 

A 

degrees of freedom- are deleted and the sign is changed on the dnf’s. 

A potential user may wish to input the eigenvectors and eigen- 
values using formatted data from their vibration package. If these 
data are in the order given above in this note, it will be necessary 
for the user to modify only the read statements for the A1 and B1 data 
blocks to reflect the desired formats. These read statements are near 
the 80th line of the FLUT subroutine of the computer program. It is 
more likely that the user’s data are in a form other than that given 
above. In this case, it will be necessary to alter either the data 
or the program. 

(b) Calculated versus interpolated aerodynamics.- 

Data Blocks B, P, and Q differentiate between the reduced fre- 
quencies used for aerodynamic calculations (RFQ; and the reduced 
frequencies used for flutter analyses (RFT). The motivation for 
these two sets is explained in Section III and has mainly to do 
with minimizing the utilization of computer resources. Typically, 
the aerodynamics are computed explicitly at three reduced 
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frequency values and these results can be used to Interpolate to 
many reduced frequencies. The interpolation has worked well on 
several examples, with the following conditions: 

(i) It is advised that extrapolation of the aerodynamics be 
kept to a minimum. The program will produce answers, 
but their quality may be suspect. 

(ii) Interpolation using k values that differ by an order of 
magnitude should always be acceptable (unless extremely 
accurate results are needed such as, for example, those 
required for flutter derivative information). Ranges 
larger than an order of magnitude should be used with 
caution, 

(iii) If an RFT equals an RFQ,no interpolation is performed. 
Likewise, if NRF <; 3, no Interpolation can be performed, 

(iv) An instance where interpolation can be extremely valuable 
is when added definition is needed for a particular 
reduced frequency range. If the aerodynamics from the 

RFQ values have been saved on tape, they can be read off the 
tape and the flutter analysis done at the added reduced 
frequencies without a new aerodynamic analysis. 

(c) Panels, strips and boxes, - 

The terminology used In the definitions of the aerodynamic 
representation is sufficiently confusing to warrant some 
clarification. Figure 3, which is another representation of the 
wing of Figure 2, is used for illustrative purposes. 


20 


A panel is a large aerodynamic surface which can be represented 
by a trapezoid. The edges of the panel must be parallel to the 
free stream direction. There are two panels in Figure 3, 

A strip is a chordwise segment of a panel. It extends from 
the leading to the trailing edge. There are six strips in Figure 3. 

Each strip can be further subdivided into boxes. In Figure 3 
each strip has four boxes. 

Reference 2 has a number of rules for setting up the aerodynamic 
model. Among these are: 

(i) The aspect ratio of the boxes should be unity or less, 

(ii) There should be a minimum of four boxes per strip, with 
even more at high reduced frequencies. 

Both of these"rules“ have been violated by the author while performing 
analyses with little adverse effect, A brief study to determine 
how refined the aerodynamic grid had to be also indicated that the 
rules were conservative, particularly at low k’s. 

The use of AIC‘s involves collocating the aerodynamics from 
each strip at a chordwise location given by data Blocks M and N 
and a spanwise location at the center of each strip. 

(d) The variable NCORE.- 

The size of the problems under study is given by NCORE. This 
parameter dictates which of two subroutines is to be used to solve 
the linear systems of equations in the aerodynamics package. If the 
value of NCORE is less than 3700, all the calculations can 
be performed in core. If it is greater than 3700, auxiliary storage 
is required and a separate subroutine handles this case. Unfortunately, 
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the equations are solved with a change in sign between the two 
subroutines. Therefore, if a case is run with NCORE greater 
than 3700, the sign of the AlC's must be changed before they are 
used in FLUT. A simple way of doing this is to input values 
of RHON in data Block T that are the negative of the actual 
density ratios. 

(e) Very low reduced frequencies.- 

If static divergence types of instability are being investigated, 
very low values of k (reduced frequency) should be used. In fact, 
static divergence only occurs, by definition, at k = 0. 
Unfortunately, the aerodynamics package is configured in a way that 
gives different types of AlC's at zero k than it does at nonzero 
values. These static AlC's are not compatible with FLUT. Therefore 
static divergence must be investigated in the limit as k ->• 0. 

The crossover from dynamic to static AlC's occurs at k = 1x10'^, 
representing a lower limit for k's in this program. A secondary 

r 

crossover occurs when k/FL (see Block D) is 5x1 o". For values of 
k less than this value, only steady state aerodynamics are 
calculated with no increment due to oscillations. This should 
not cause any problems, but users should be aware of this transition 
when they are Interpolating across this value, 

VI, Program Execution 

If the input described in the previous section Is correctly entered, 
results will be printed out at the values of reduced frequencies and 
atmospheric densities specified. It is the intention of this section 
to explain briefly the algorithm used to achieve these results and to then 
provide some information that could aid in the Interpretation of the results. 
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(a) Algorithm for FLUT 


Figure 4 sfiov/s a simple block diagram for the algorithm. 

The diagram indicates that the program is initiated by reading 
the basic parameters of Block B. The broken line to the box 
"calculate AIC's" indicates that this calculation, which uses input 
from Blocks C through 0, is performed only if the results have not 
been saved on tape during a previous run. After the remaining input 
is read from cards, the vibration and aerodynamic data is read from 
tape and the generalized aerodynamic forces at the RFQ reduced 
frequencies are calculated. The complex eigenvalue analysis is then 
performed for all the RFT reduced frequencies and RHON density values, 
with density on the inner loop and reduced frequency on the outer loop. 

The program is contained in two overlays plus a zero overlay that 
controls the program execution. 

Overlay one contains the flutter analysis ahd has f ve subroutines 
of its own. Four of these are a "black box" that perfor is the 
complex eigenvalue analysis. 

Overlay two contains the unsteady aerodynamics analysis. With 
slight modification, this is the program developed by Giesing, Kalman 
and Rodden (Ref, 4 ) . 

(b) Interpretation of the Results 

The aeroelastic stability analysis contained in this program 
Is based on the V-g method of flutter analysis. Any standard 
reference text (e.g. Ref. 1) on aeroelasticity contains a description 
of this method. The principle features of this analysis are shown 


in the sketch of Figure 5 which depicts two branches of a V-g curve at 
a constant atmospheric density. FLUT solves the complex eigenvalue 
problem at each reduced frequency and density. It is the user's task 
to take the output, which gives values of velocity, V, artificial 
damping, g, and frequency and plot curves like the one shown in 
Figure 5. The plotting of these curves is something of an art with 
some helpful tips listed here: 

(1) There is a branch on the V-g curve for each of the normal modes 
included in the stability analysis. Typically, only a handful 
are in the velocity range of interest with the remainder at 
physically unrealistic velocities. 

(2) For a given branch, high k values correspond to low velocities and 

vice versa: i.e., in Figure 5, k- < k. There can be exceptions 

1 1+ j 

to this where a branch doubles back on itself slightly. 


(3) The frequency results serve as an aid in tracking a branch through 
a range of k values. At very high values of k,the velocities 
approach zero and the frequencies approach the natural vibration 
frequencies. Frequently, it is possible to start from a high k 
value, identify the branches and follow them by assigning results 
at lower k's to branches that have a corresponding frequency. 

(4) When in doubt, , results at intermediate k values should be obtained. 
If the aerodynamics at the'RFQ values have been saved on tape, it 1s 
simple and inexpensive to run more RFT^s.^ With enough data points, 
the correct branch assignments should be obvious. Once the branches 
have been drawn, the only remaining task is an assessment of the 
stability of the system. There are two types of instability that 

the wing can exhibit. The first, type is depicted 
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in Figure 5 by one of the branches crossing the V-g axis. In this 

case, at g = 0 the system is said to be in neutral equilibrium 

at the corresponding velocity. Given a disturbance, the system 

will oscillate harmonically at the frequency that corresponds to this 

point. At higher velocities g is positive and the system is 

unstable with a response that is characterized by divergent oscillations. 

A second type of instability that sometimes occurs is manifested by 
a V-g branch that approaches g = 0 from below, but never crosses into 
positive g. In the limiting case as k + 0, this phenomenon corresponds 
to static aeroelastic divergence. The response of the wing in this case 
is characterized by a nonoscillatory divergence. 

As mentioned, the V-g curves are drawn at a constant atmospheric 
density. The effect of altitude on the flutter velocity can be studied by 
plotting V-g curves at a series of RHON's. 

FLUT is programmed so that all results of a given run are obtained 
at a constant Mach number. This means that, strictly speaking, the 
results are valid only at one velocity on the V-g curve. It is extremely 
unlikely that this velocity will correspond to the speed of instability 
so that the results are inaccurate. If greater precision is desired, it 
is necessary to run the program at several Mach numbers and then interpolate 
to find the point where Mach number and the flutter or divergence velocity 
correspond. An example of this is shown in Figure 6, which has been taken 
from Reference 9 . In this figure, flutter velocities calculated at three 
altitudes and Mach numbers are given by the square data points. At a 
given altitude, the flutter point is seen to decrease slightly with Mach 
number. The three lines emanating from the origin are simply lines where 
the velocity and Mach number agree for each of the altitudes. The circled 
data points are therefore matched flutter points where the flutter 
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velocity and Mach number agree for the given altitude. It should be 
noted that, in this case, the matched flutter point for the two highest 
altitudes are at Mach numbers that are outside the range where the 
extrapolation from calculated results is quantitatively meaningful. 

VII. Example: Go! and’ s wing 

Reference 10 is a journal article that presents an early technique used 
for flutter analysis. In the article, the author, Martin Goland , presents 
a simple example that can be used to verify the FLUT program. Figure 7 
depicts the cantilevered wing he used in his example and lists the relevant 
structural parameters. 

A vibration analysis of this model was performed using the PASS computer 
program (Ref. 9). The 11 structural grid points used in the analysis were 
arranged along the elastic axis at the following intervals: (given as a 

fraction of the semi-span): 0.0, .0625, .1875, .3125, .4375, .5625, .6825, 

.8125, .9167, .979, 1.00. These points, except for the first and the last, 
are at the center of the aerodynamic strips selected for the unsteady aerodynamic 
analysis. It should be noted that the vibration analysis was done in inches 
so that the mode shapes are in the units of inches and radians. 

The input deck used to obtain the output listed for this example is 
listed below. The Input is identified by reference to the data blocks listed 
in Section IV.. In order to limit output, only one reduced frequency value 
was used. 

(a) Input 
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(b) Output 

The output obtained for a single reduced frequency for this example 
is given here. For the sake of brevity, some of the less meaningful output 
has been deleted. 

Figure 8 shows th4 complete V-g diagram obtained using this model. The 
figure indicates a flutter type of instability at 450 ft/sec. and a divergence 
instability at 1000 ft/sec. 
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REROD VN0niC5 FQE GOLRND WING 


flRRfiV OF REDUCED FREQUENCIES 


. 400000 

REF. CHORD = F2. 00000 REF. SEMI-SPRN = £40.00000 REF, fiREPi = 34560.00000 

M0CH NO. = 0.00000 BETR = 1.00000 


iPKJK PRNEL NO. 1 INPUT t^ALUES 


XI = 

0 .;000000 

X2 = 

72.000000 

Y1 = 0.000000 

Zl = 

0.000000 

X5 = 

0 - 000000 

K4 = 

72.000000 

V2 = 248.000000 

Z2 = 

0 .000008 

NC = 

3 NS = 10 

NDELT 

= 1 HO. 

OF PANELS = 1 




3 CHORDWISE DIVISIONS FOR PANEL 1 


^ Ci 
tr. 0 

se 

® o 

rO 

s ^ 

S’® 


0. 

.500000O0E+00 


. 10000000E+01 





10 SPRNWISE 

DIYI 

3 1 OHS FOE PANEL 

1 



0. 

. 12580000E+00 


.25000003E+00 

.37500000E+00 

. 50000000E+00 

.02500800E+03 

.75000000E+00 

.87500000E+00 


. 95833300E+CI0 

. 1000u000E-f-01 
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XT ELEMENTS 

FOR PANEL I’KY’ic 



0. 

.36000000E+02 


. 72000080E+02 

0. 

.360D0000E+02 

. 72000000E-1-02 

0. 

.36000O00E+02 


.72000000E+02 

0. 

.36000000E+02 

.720G003OE+f^£ 

0. 

. 36000000E+02 


.72008000E+02 

0. 

.3G000000E+02 

. 7£0800'30E+O2 

0. 

.36000000E+02 


.r£000000E+02 

0. 

.36000000E+02 

.72000000E^«2 

0. 

.3G000000E+02 


.7200O000E+02 

0. 

.36000000E+rj2 

.72O0BO0U£+O:=' 



to 

ETA ELEMENTS 

FOR PANEL MK-ioc. 



0. 

0. 


0. 

.300O80O0E4-02 

.30000000E-I-02 

, 30000000E+02 

.60000008E+02 

.60O00000E+02 


.60000000E+02 

. 3000000OE-[-02 

. 9OQ30000E+02 

.OBO00000E+02 

. 12000000E+03 

. 12000000E+03 


. 12000000E+B3 

. 1530800OE+03 

. 15000000E+03 

. 1500U000E-1O3 

. 1800000QE+03 

. IS00O000E+03 


. 1S000000E+03 

.21000000E+O3 

.21000000E+03 

.21006000E+33 

.2239S09EE+83 

.22999992E+03 


. 22999992E+03 

.24000O0OE+03 

. 24000000E+03 

. 24000BOOE+03 




30 ZEE ELEMENTS 

FOR PANEL HO. 1 




additional output omitted 


STABILITY DEEn/'ftTIl''E3 FOR THE. 9 STRIPS 


STRIP 

VOS 

PLUNGIN6 

NORMAL 
SPAN LOAD 

DEEIkATU^ES 

NORMAL 
SPAN MOMENT 

NORMAL 

AILERON HINGE 

MOMENT 

NOEMAI, 

TAB HINGE MOMENT 

I 

.0625 

-3.84609 

-.41768 

.3HC95 

17183 

0,00000 

0.00000 

0.00000 


0.060O( 

2 

. 1375 

-3.83344 

-.42330 

,32064 

17091 

0.00000 

0.00000 

0.00600 


B.0000( 

3 

.3125 

-3,30285 

-.45162 

.31936 

16897 

0.00000 

0.00000 

0.00000 


3. 00001 

4 

.4375 

-3 , 74522 

-.43699 

.31828 

16576 

0.00000 

0.00000 

6.00000 


3 . 00001 

5 

.5625 

-3.641 17 

-.53482 

.31519 

16087 

0.00000 

0.00000 

6.00000 


0 . 00001 

6 

.6375 

-3.45222 

-.59094 

.30900 

15355 

0.00000 

0,00000 

0 , 00000 


0.60B0I 

7 

.8125 

-3.09784 

- . 63654 

. 29520 

1419G 

8.00000 

0.00000 

0 . 00000 


0.00001 

B 

.9167 

-2.46180 

-.60351 

. 25906 

12211 

0 . 00000 

0 . 00000 

0.80000 


8.00001 

9 

.9792 

-1.57843 

-.43311 

. 18039 

08393 

0 . 00000 

0 . 00000 

0.00000 


0 .00001 


C-SUB-L OR V 
C-SUB-tl OR M 
C-5UB-M OR H-DEL-fi 
C“SUM-ri OR H-DEL-T 
C-SUB-SOEIPT-L 


"3.4483<i>S3S 
-.83336067 
0 . 00000000 
0 . 00000000 
0 . 00000000 


-.51176213 
-.32428550 
0.00000000 
0 . 00008000 
0.00000000 


0 . 00000000 
0,00000000 
0 . 00000000 
0.00000000 


0,00000000 
0.00000000 
0 . 00000000 
0.00000000 


PITCHING DERI1-^1=)TI1/ES 


STS IP 

Y05 

MOEMPIL 
SPAN LORD 

HOSMRL 
SPRH MOMENT 

NOKMRL 

flILEEON HINGE 

MOMENT 

NOEMOL 

TRB HINGE MC-: 

•TENT 

1 

.0625 

-3.87086 

-1.72064 

.37339 

-.25285 

0 . 00000 

0 . 00000 

0.00000 

O.0600E 

2 

. 1075 

-3.B53S7 

-1.72792 

.37275 

-,29203 

0,00000 

0 . 00000 

0.00000 

0.O0O0E 

3 

.3125 

-3.S155S 

-1.74184 

.37141 

-.25024 

0 . 00000 

0.00080 

0.00088 

0.OG50C 

4 

.4375 

-3.74576 

-1.76032 

.36884 

-.28734 

0 . 08000 

0.00600 

0.00600 

0 . 000.0E 

5 

.5625 

-3.62478 

-1.77746 

■ .36426 

-.28258 

0.00000 

0.00000 

0.00060 

0.00087 

6 

.6875 

-3.41457 

-1.77733 

.35592 

-.27628 

0.00000 

0 .08000 

0.00000 

0. 8000E 

? 

.8125 

-3.03805 

-1.71551 

. 33884 

-.26406 

8.00000 

8.00000 

0 . O000B 

0 , 08300 

8 

.5167 

-2.35365 

-1.48026 

,29703 

-.23435 

0.00000 

0.00000 

0,00000 

0.0O00E 

3 

.5752 

-1.52555 

-1.00642 

,20680 

-. 16444 

0.00000 

B. 00000 

0 . 08080 

0.080BC 


CjO 

CD 


C-SUB-L OR Y 

-3 . 43348492 

-1 .69291530 

0 . 00000000 

0 .00000008 

C-5UB-M OR N 

-.T«149973 

-.83320968 

0.00000000 

0.00000000 

C-5UB-M OE M-DEL-FI 

0.08000000 

8.00000000 

0.00000000 

0.00000000 

C-SUM-n OE N-DEL-T 

0.00000800 

0. 00000008 

O.00G00000 

0,08000000 

C-SUB-SCRIPT-L 

Q. 00000000 

0 , 00000000 




PFES5UEE COEFFICIBNtS FOE MODE NO, I 


STRIP HO. 

xoc 

C-SUB-P/ 


1 -r 

.1250 ■ 

“.00216055 

-.06469893 

' ■ - 1- 

.6250 

.01144241 

-.02073304 

S : 

.1250 

-.00197410 

-.06451118 


.6250 

.01150291 

-.02087634 

■ 3 :■ 

.1250 

-.00156339 

-.06407571 

3 ■ 

.6250 

.01162444 

. -. 02043216 

: 4 

, 1-250 , 

-.60098209 

-.06324380 

4 

.6250 

.01180462 

-.01997735 

5 ;■ 

. 1250 

-.00013784 

-.06 174005 

&/ 

.6250 

.01202273 

-.81916688 

i5 

. 1250 

.00092340 

« 

-.05399584 


.6250 

.01220850 

- .01772024 


, 1250 

,00203641 

-.05373628 

7 

.6250 

.01210397 

-.01510458 

8 

.1253 

,00248564 

-.04379054 

8 

. 6250 

.01092577 

-.01091605 

St: - 

. 1250 

.00134846 

-.028 ,’1232 

■ a 

. 6250 

.00767611 

-.00636440 


;p NO. 

F RESSURE COEFF 1C I EMTS 
XOC 

FOR MODE MO. £ 

C-3UB-P 

1 

. 1250 

-6.061 IS 109 

1 

.6250 

-1,68053009 

o 

. 1250 

-6.03872976 

2 

.6250 

-1,66901027 

3 

. 1250 

-5 . 98799639 

3 

.6258 

-1,64310504 

4 

, 1250 

-5 . 89536097 

4 

.6250 

-1.59615691 

5 

. 1,250 

-5.73429186 

5 

. 62.50 

-1.51526464 

6 

• 1250 

-5 , 45205670 

6 

.6250 

- 1 . 376PSC 43 

7 

. 1250 

-4.94024193 

7 

.6250 

-1. 135SG651 

8 

, 1250 

-4 01267299 

8 

,6250 

-.77470418 

9 

. 1.250 

-2.62731055 

9 

.6250 

-.42409557 


EIGENI^LUES 

3322. 

9163. 

. 6034E+05 
. 122SE+S6 
. £08?E+0€ 

. 4'364E+0S 
. r45SE+0G 
. 9209E+06 
. 1214E+07 
. 1762E-1-07 

EIGEHI-^CTOKS 


DOF 1 2 

1 -.I125E-0I .327i3E-0Z 

2 -.I477E-02 ' .69S3E-03 

3 .884 IE-03 . ;061E-P2 

4 -.95PQE-01 .4598E-8I 

5 ".4P31E-R2 


.3227E-01 
.405SE-02 
[3S8E-91' 
. 2254 
.7rO”2-n2 


4 

.5578E-01 
.7077E-02 
.P027E-0 ■ 
,40'^ 3 
. 14b2E-01 








-.8588 1988 
-2.53246845 
-.87081777 
-2.58501337 
-.89440152 
-2.58927463 
- . 9278 1627 
-2,53282148 
- . 96548933 
-2.53943133 
- . 99213325 
-2.56251550 
-.96804369 
—2 . 46297629 
-.80931389 
-2. 151 19798 
-.52981680 
-1.4S30226S 


5 


€ 


.4S78E-01 -.798IE-02 .204RE-C1 . :9j2 

,5981E-riE -.1165E-0P .2189E-82 . 2"05 

-.2645E-0J -.47e3E-01 -.9'-»409-81 . 

• 3054 — .3773E-01 i2i:,4E-i:l I.G.- ■ 

.O223E--0: -.792'i'" 1.2 .:’;24i-L)3 . 2 ■• 7. F- “ 
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... y.,,: 

-'►2476 

.1236 ^ 

.4557 

.8439 


-'.eS36E-02 

.3150EI-02 

•6739E-02 

.1330E-01 

'• 9 ■ 

►433BE~8Z 

.3327E-01 

-.48S7E-01 

.5777E-0! 

10 

: -.4535 

.2331 

.5955 

1.124 

,11 

-,7572E-^02 

' .4M4E-02 

.2184E-02 

, 4329E-02 

12 

.5929E-0Z 

■ . 4454E-81 ■ 

-.43I2E-01 

..43i7E‘-0I 

13 

-.6974 

" .3632 r 

.5737 

1.059 - 

14 

>r.8616E-02 

; .4852E-82 

- . 3677E-02 

-. 9052E-B2 

19 

.7348E-02 

.5395E-01 

-.2180E-01 

. M71E-01 

16 

-,9663 

.5217 

.3821 . 

.5329 

17 

-.9246E-02 

.S343E-82 


-.2223E-01 

18 

.852 IE-02 

.6120E-01 

.9029E-02 

- . 1227E-0 1 

19 

-1.249 

.6864 

.6375E-01 

-.2346 

zd 

-.9547E^2 

.5600E-02 

I201E-O1 

-.3126E-01 

21 

.9375E-02 

.6608E-O1 

.3930E-0I 

-.2332E-01 

22 

-1.439 

.8276 

-.2514 ■ . 

-1.862 

23 

-.9S28E-02 

.5674E-02 

-. 1239E-01 

-.3429E-01 

24 

.979BE-^02 

.6827E^01 

.5633E-01 

-.2410E-01 

25 

-1.633 

.9127 

-.4474 

-1,580 

26 

-.9635E-0E 

- 5681E-02 

-. 1310E-01 

-.3463E-01 

27 

.9B9.iE-02 

.6877E-01 

.6068E-0! 

-.2340E-31 
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VIII, Concluding Remarks 


The preceding has pointed out the main features of the FLUT program. It 
has been successfully applied to several example problems and to the analysis 
of the Lockheed Subsonic Oblique Wing Transport Concept (Raf. 9). Recommendations 
for enhancements and conclusions reached during its use can therefore be made 
here. Experience with the use of this program has shown that, in performing a 
flutter analysis, it is extremely helpful to use two separate programs to do the 
same analysis. While it may seem like a luxury to do the same thing twice, in 
practice one program serves as a check on the other, thereby pointing out 
errors in input and finally adding confidence to the results when matched 
results are obtained from the two programs. For instance, while using FLUT, 
it was found useful to perform the same analysis using NASTRAN (Ref. u). 

A minor recommendation is that the unsteady aerodynamics package is several 
years old and could probably be improved upon. There may, in fact, be newer 
packages that are more suitable. Alternatively, it may be v/orthwhile to per- 
form some re- programming of the algorithm. The package was designed for a 
machine with limited core space and makes liberal use of tapes for intermediate 
storage. This results in awkward data handling that could certainly be stream- 
lined. The recommendation is minor in the sense that it is not clear that the 
benefits would justify the expense. FLUT has already been developed so 

that it minimizes the calculations in this most inefficient part. 

A more ambitious recommendation is that this flutter analysis program has 
many features which make it attractive for use as a design tool. With the 
major addition of flutter derivative calculations, the program could be used 
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in optimization studies with constraints on the aeroelastic stability. 

References 12 and 13 present various ways that the program could be used to 
achieve this. 

Conclusions regarding the use of aerodynamic interpolation and the reduction 
of rigid body modes as well as other concepts are made throughout the text. 
Aerodynamic influence coefficients are felt to be of sufficient merit to 
reiterate their usefulness a final time. Finally, the program has been shown 
to be efficient and reliable in its execution. 
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Figure 7.- Goland's Uniform Cantilever Wing. 
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